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AIM 


TtiO  terrain  analysis  system  of  programs  derives,  from  a  matrix  of  altitude 
data,  the  descriptive  characteristics  of  gradient,  aspect,  and  profile  and  plan 
convexity  for  the  central  point  of  each  3x3  square  set  of  adjacent  grid  points. 
The  characteristics  are  derived  algebraically  from  a  quadratic  fitted  by  the 
least  squares  method  to  each  set  of  9  points.  For  the  whole  matrix,  the  mean, 
standard  deviation,  skewness  and  kurtosis  are  found  for  the  calculated  altitude, 
gradient  and  profile  and  plan  convexities.  The  vector  mean  of  the  aspect  is 
found  and  also  the  regression  of  each  of  the  ocher  characteristics  with  the  sine 
and  cosine  of  the  aspect.  The  descriptive  characteristics  are  displayed  in 
several  ways,  viz,  histograms,  scatter  diagrams  and  line  printer  density  maps. 

A  map  is  produced  on  the  graph  plotter  giving  the  direction  and  magnitude  of 
the  gradient  at  each  grid  point.  The  programs  are  written  in  FORTRAN;  they 
call  SPSS  and  MIDAS  routines  and  run  under  the  MTS  system  on  the  NUMAC  IBM 


370/168  computer.  |k  The  system  and  the  special  programs  were  written  by  Margaret 
and  Tony  Young,  to  Specifications  by  Ian  S.  Evans.  Reasons  for  these  specifi¬ 
cations  are  given  in\the  Final  Report. 

MATHEMATICAL  DERIVATIONS 


Input  data  The  program  has  been  designed  to  accept  a  rectangular  matrix  of 
altitude  data.  It  assumes  that  the  data  are  derived  from  points  equally  spaced 
in  the  north-south  and  east-west  directions,  that  the  first  row  is  the  northern¬ 
most,  and  that  it  is  to  be  read  from  west  to  east.  The  program  requires  the 
distance  in  metres  between  grid  points,  the  number  of  columns  in  the  matrix  of 
data  and  a  scaling  factor  to  convert  the  altitudes  to  metres.  For  an  irregular 
area  such  as  a  drainage  basin,  the  remaining  grid  points  within  an  enclosing 
rectangle  are  given  as  zero  and  any  set  of  3x3  points  which  contains  a  zero  is 
ignored , 


Least  squares  equation  For  each  3x3  square  of  grid  points  a  quadratic  equation 
is  fitted  by  the  method  of  least  sum  of  squared  deviations.  Let  this  be: 

z  =  ax2  +  by2  +  cxy  +  dx  +  ey  +  f 


and  let  the  mesh  points  of  the  grid  be  g  metres  apart.  Then  the  least  squares 
solution  is  obtained  by  solving  the  standard  matrix  equation: 


T  T 
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where  a  is  the  6x1  matrix 


and  z  is  the  9x1  matrix  of  heights  at  the  9  grid  points  (~g,g) 
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The  matrix  F  is  thus  given  by: 
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As  the  matrix  is  only  evaluated  once  it  is  neater  to  use  this  form  rather 
than  to  compute  it  and  store  the  coefficients. 
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Derivation  of  descriptive  characteristics  The  gradient,  aspect,  profile 
(or  vertical)  convexity  and  plan  (or  horizontal)  convexity  are  derived  from 
the  quadratic: 

z  =  ax2  +  by2  +  cxy  +  dx  -rey  +£ 

Also  if  x  =  rcosG)  ,  .  , 

.  this  transforms  to 
y  =  rsir.Q) 

z  =  r2(acos“©  +  bsin20  +  csinQcosG)  v  r(dcos0  +  csinG)  +  f 
The  gradient  of  the  origin  (central  point)  is  given  by: 
dz 

(t— )  A  =  2r(acos2G  +  bsin20  *  csinGcosG)  +  dcosG  +  esinG 
or  r=0 

=  dcosG  +  esinG  (=  G  say) 

•  qG  •  •  £ 

For  max.  or  min.  -rr  =  0,  i„e. -dsinQ  +  ecosG  =  0  ,or  tanG  =  •? 

dG  *  d 

which  gives  the  aspect  angle,  G0 

The  standard  two-dimensional  expression  for  curvature  is: 


>  < 


Using  the  convention  that  convex  surfaces  have  positive,  and  concave  surfaces 
have  negative  curvature,  this  gives: 


Plan  curvature  (plan  c) 


and  Profile  curvature  (prof  c) 


-d2r 

dz2 


for  prof  c 


-  k  - 


which  gives: 

prof  c  = 
‘ (r=0) 


-2(acos20  +  bsin20  esinOcosO) 
{1  +  (dcosO  +  esin0)2}  '2 
-2(ad2  +  be2  +  ced)/(e2  *  d2) 

{ 1  +  (d2  +  e2)  / (e2  +  d2) }  ^2 

—2  (ad2  +  be2  *  ced) _ 

(e2  +  d2)(l  +  d2  +  e2)  ^2 


For  plan  c 


ax2  +  by2  +  cxy  +  dx  +  ey  +  f  =  const. 


differentiating  with  respect  to  x  gives: 

2ax  +  2by~  +  cy  +  cx~-  +  d  +  e^—  =  0 
J  dx  J  dx  ax 

.  dy  _  -<2ax  +  cy  +d) 

2by  +  cx  +  e 


At  (0,0)  &  «.  ^ 
*  dx  e 


d2y  =  L{ 

dx2 


(2ax  +  cy  +  d)(2b^~  +  c)  -  (2by  +  cx  +  e)(2a  +  ■—£) 
(2by  +  cx  +  e)2 
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Now  plan  c 


*  (ft)2/2 
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For  a  function  z=f(x,y)  to  have  zero  gradient 


3z  A  ^  .  3z  A 

_,°  ra  _  =  o 

For  the  quadratic  ax2  +  by2  +  cxy  +  dx  +  ey  +  f  =  z,  zero  gradient  occurs 
at  a  point  x,y  given  by: 

3z 


3x 


oZ 

3y 


=  0  =  2ax  +  cy  +  d 


=  0  =  2by  +  cx  +  e 


i.e. 


ec  -  2bd 


y  = 


dc  -  2ae 


4ab  -  c2  4ab  -  c2 

The  condition  for  a  saddle  point  is: 

<  o 

3x2  3y2  3x3y 

i.e.  2a  .  2b  -  c2  <  0 
i.e.  4ab  <  c2 

For  a  maximum  or  a  minimum: 
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gradient  points  awdj-  from  the  point  (x-0,y=O)  to  which  all  the  calculations 
in  the  analysis  program  refer. 


Description  of  zero  gradient  points  At  the  central  point  in  some  cases 
.(most  commonly  for  lakes)  the  values  calculated  for  d  and  e  are  both  zero, 
i.e.  the  gradient  is  zero  and  the  aspect  is  vertical,  and  thus  the  aspect 
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angle  is  indeterminate.  As  a  complete  set  of  characteristics  does  not  exist 
for  such  a  point,  no  values  are  stored  for  inclusion  in  the  main  statistics. 
It  is  possible  to  find  the  shape  of  terrain  at  such  a  point  by  calculating 
profile  convexity  and  the  maximum  and  minimum  values  determine  this  shape  as 
follows  1 


This  can  be  written:  2(Acos(20  +  B)  +  C) 

where  Acos(20  +  B)  +  C  =  A(cos20cosB  -  sin23sinB)  +  C 

=  AcosB(cos20  -  sin20)  -  AsinB(2sin9cos0)  +  C(sin20  +  cos20) 
=  A(cosB  +  C)cos20  +  (C-AcosB)sin2©  -  2Asin3cos0sin0 


i.e. 


giving 


a  -  AcosB  +  C 

b  =  C  -  AcosB 

c  -  -2AsinB 

„  _  a  +  b 
C  2 

AcosB  =  --y- 
AsinB  =  -  ~ 


giving 


A2  =  (• 


a  -  b.2  ju  (S.\ 2  -  (a  ~  b)~ 

2  }  ‘  K2J  4 


and  tan  B  « 


For  gradient  =  0 
nrof  c  =  - 

dr2 

=  -2(Acos (20  +  B)  +  C) 

Max  prof  c  =  -2C  +  2 1 A. J  =  -(a  +  b)  +  ((a  -  b)2  +  c2)2 
Min  prof  6  =  -2C  -  2  j  A.}  ~  -(a  +  b)  -  ((a  -  b)2  +  c2)2 

For  the  cases  when  the  surface  is  a  horizontal  plane,  i.e.  a  plain  or  lake, 
a=b=c=d=e=0  i.e.  A  =  C  =  0 
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For  a  sunrp.it,  i.s.  convex  surface 
prof  c  >  0  for  all  0 

i.e.  | A j  <  |cj  and  C  <  0  c2  <  4ab  and  a  +  b  <  0 

For  a  ridge 

prof  c  £  0  for  all  0 

i.e.  J  A |  =  | C |  and  C  <  0  c2  =  4ab  and  a  +  b  <  0 

For  a  saddle 

prof  c  has  positive  and  negative  values 
i o  e .  | A  j  >  [ C |  c2  >  4ab 

For  a  valley 

prof  c  £  0  for  sil  0 

i.e*  | A |  =  J C j  and  C  >  0  c2  =  4ab  and  a  +  b  >  0 

For  a  pit 

prof  c  <  0  for  all  0 

i.e.  J A |  <  jcj  and  C  >  0  c2  <  4ab  and  a  +  b  >  0 

Thus  the  maximum  and  minimum  values  of  prof  c  give  the  shape  of  the 
terrain  as  follows: 


max  >0 

min-  >0 

summit 

max  >0 

min  =0 

ridge 

max  >0 

min  <0 

saddle 

max  =0 

min  =0 

plain 

max  =0 

min  <0 

valley 

max  <0 

min  <0 

p:  : 

Other  methods  of  detecting  topographic  features  are  given  by  Peucker  and 
Douglas  (1975)  and  Grender  (1976), 


wtsszxny:;- 
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PROGRAM  DESCRIPTIONS 

The  main  program  The  main  program  fits  the  least  squares  quadratic  and 
calculates  the  descriptive  characteristics  gradient,  aspect  and  the  profile 
and  plan  convexities.  Initially  the  program  reads  in  and  scales  two  rows 
of  altitude  data,  and  subsequently  one  row  at  a  time.  Using  three  rows  it 
tests  for  zero  values  then  calculates  the  characteristics  for  each  point  in 
the  row.  If  a  test  for  zero  gradient  is  positive,  the  altitude  and  maximum 
and  minimum  values  of  profile  convexity  are  printed,  with  a  description  of 
the  type  of  point.  Such  points  are  not  included  in  the  main  statistics  as 
values  of  the  characteristics  are  not  complete.  However,  sums  of  squares 
and  higher  products  of  altitude  are  stored  separately  for  additional 
statistics. 

The  characteristics  for  non-zero  gradient  points  are  stored  in  a  file 
for  use  by  other  programs.  The  sum  of  squares  and  cross  products  are 
calculated  for  the  statistics »  Double  length  working  is  used  for  the  sums 
of  squares  and  higher  products  as  loss  of  accuracy  is  possible  when  forming 
sums  of  up  to  40,000  points,  particularly  for  the  higher  moments.  When  the 
end  of  the  altitude  data  is  reached  the  program  calculates  (for  non-zero 
gradient  points)  the  mean,  standard  deviation  and  the  skew  and  kurtosis,  and 
also  correlation  coefficients  for  calculated  height,  gradient  (which  is 
expressed  in  degrees)  and  profile  and  plan  convexities  (which  are  in  degrees 
per  100m),  Mean,  standard  deviation,  skewness  and  kurtosis  are  calculated 
for  all  points  for  the  characteristics  calculated  height  and  gradient,  and 
for  plain  and  non-zero  gradient  points  for  prof  c  and  plan  c.  The  vector 
mean  of  the  aspect, and  the  vector  mean  of  the  aspect  weighted  by  the  gradient, 
are  also  calculated.  The  statistics  are  printed  and  the  means,  standard 
deviations,  maxima  and  minima  are  stored  for  use  by  other  programs. 

The  main  program  calls  several  specially-written  programs  and  parts  of 
the  SPSS  and  MIDAS  packages. 

The  mapping  programs  One  of  the  ways  in  which  the  characteristics  are 
displayed  is  a  series  of  line  printer  maps.  The  value  of  a  characteristic 
at  each  grid  point  is  indicated  by  the  density  of  printing.  A  blank  space 
indicates  a  missing  value  at  that  point,  due  either  to  the  shape  of  the  area 
or  to  a  point  of  zero  gradient  for  which  some  characteristics  do  not  exist. 

In  aspect  data  8  density  values  are  used,  the  darkest  indicating  north  and 
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the  ligate.: t  south,  while  east  and  west  directions  have  printer  characters 
different  in  form  but  similar  in  density. 

Multiple  regression  and  histograms  The  SPSS  and  MIDAS  programs  available 
in  the  NUMAC  system  are  used  for  multiple  regression,  scatter  diagrams  and 
histograms.  SPSS  is  called  to  produce  the  following  regressions: 

Calculated  height,  gradient,  profile  convexity  and  plan  convexity 

each  as  functions  of  sine  and  cosine  of  aspect. 

Gradient  as  a  function  of  calculated  height,  square  of  calculated 

height,  sine  and  cosine  of  aspect. 

Profile  and  plan  convexity  each  as  functions  of  calculated  height, 

gradient,  sine  and  cosine  of  aspect. 

These  combinations  were  chosen  as  interpretable  sets,  of  controlling 
variables,  influencing  the  (first-named)  dependent  variable. 

SPSS  is  also  used  to  produce  scatter  diagrams  of  each  pair  of  the 
following:  calculated  height,  aspect,  gradient,  profile  convexity  and  plan 
convexity. 

The  SPSS  scattergram  program  prints  the  number  cf  values  at  each  position 
up  to  a  maximum  of  9  -  some  detail  is  lost  in  this  way  but  the  shapes  of  the 
distributions  can  easily  be  seen. 

The  MIDAS  programs  are  called  to  produce  a  histogram  of  each  characteristic. 
The  limits  and  scales  to  be  used  by  this  program  are  calculated  by  the  limits 
program. 

Slope  plotting  program  The  slope  plotting  program  produces  a  map  showing 
the  direction  and  magnitude  of  the  gradient  at  each  grid  point.  For  each 
point  for  which  the  gradient  is  not  zero  an  arrow  is  plotted,  with  the  centre 
of  its  shaft  at  the  grid  point.  The  direction  of  the  arrow  indicates  the 
line  of  maximum  gradient,  and  the  shaft  is  one  of  5  different  lengths  depending 
upon  the  magnitude  of  the  gradient. 

Class'  limits  As  the  system  of  programs  is  designed  to  analyse  data  from  a 
wide  variety  of  different  types  of  terrain,  it  is  impossible  to  determine  in 
advance  the  class  boundaries  and  limits  which  will  produce  acceptable  maps 
and  histograms.  This  is  especially  true  for  altitude  where  the  variation  can 
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be  as  little  as  20  metres  or  over  1,000.  For  the  histograms  produced  by  the 
MIDAS  program  the  height  interval  is  chosen  to  produce  between  60  and  i 20 
classes,  with  a  range  from  just  below  the  minimum  to  just  above  the  maximum. 

In  the  density  maps  (apart  from  the  aspect  map  which  has  8  fixed  limits) ,  6 
densities  are  used  and  it  is  necessary  to  choose  limits  which  give  similar 
numbers  of  points  for  each  density  value.  For  altitude  and  gradient  the 
points  of  division  between  classes  use  the  mean  and  standard  deviation  (  sd)  as 
follows : 

mean  -1.2sd  ,  mean  -0.6sd  ,  mean,  mean  +0.6sd  ,  mean  +1.2sd 

which  for  a  normal  distribution  give  a  fairly  even  distribution  between  classes 
(Evans,  1977). 

For  gradient,  a  map  using  fixed  limits  was  also  produced  to  facilitate 
comparison  of  different  matrices. 

Organisation  The  complete  analysis  of  an  altitude  matrix  requires  several 
programs  to  be  called  using  data  mainly  produced  by  the  main  program.  The 
SPSS  and  MIDAS  programs  require  a  specified  format  for  the  data.  This  was 
achieved  by  creating  a  master  file  containing  all  the  run  commands  and  all 
the  information  for  each  program.  An  organisation  program  to  calculate  the 
limits  and  class  boundaries  was  written.  This  program  reads  data  from  a 
specified  file,  calculates  the  limits  and  writes  these  back  into  the  master 
file  in  the  required  format.  As  several  programs  require  the  same  data,  i.e. 
title,  file  names,  etc.,  these  are  read  into  part  of  the  master  file  and 
picked  up  as  required.  These  are: 

Title  of  run 

Name  of  file  containing  altitude  matrix 
Name  of  file  for  detailed  results  -  this  is  usually  dummy 
Name  of  file  (possibly  temporary)  for  storage  of  characteristics, 
gradient  aspect,  etc,,  to  be  used  by  other  programs 
Name  of  file  for  limits 
Grid  size  in  metres 
Scaling  factor  for  heights 
Number  of  columns  in  altitude  matrix 
Format  of  data 

The  files  for  results  and  limits  must  be  empty  initially;  if  the  same 
files  are  used  for  several  runs  they  must  be  emptied  before  each  run. 


PROGRAM  DETAILS  AND  RUNNING  INSTRUCTIONS 


Main  program  OMY8  The  main  program  accepts  a  rectangular  matrix  of  altitude 
data  together  with  a  scaling  factor.  Using  the  standard  least  squares  method 
it  fits  a  quadratic  to  sets  of  3x3  adjacent  points,  and  from  the  quadratic  it 
derives  the  calculated  altitude,  gradient,  aspect,  profile  convexity  and  plan 
convexity  for  each  point.  For  the  whole  matrix  the  mean,  standard  deviation, 
skewness  and  kurtosis  of  the  characteristics  are  calculated,  as  are  correlations 
between  calculated  altitude,  gradient,  profile  convexity  and  plan  convexity. 

Input:  Cards  A.- 1 
A  Data  header  (20A4) 

B  Name  of  file  containing  altitude  matrix 

C  Name  of  file  for  output  of  detailed  results  (usually  DUMMY) 

D  Name  of  file  for  storage  of  calculated  values  for  use  by  other  programs 

E  Name  of  file  for  storage  of  maximum  and  minimum  values  for  use  by  the 

limits  program 

F  Distance  in  metres  between  grid  points  of  altitude  matrix 
G  Scale  factor.  Each  value  in  the  altitude  matrix  is  multiplied  by 
this  scale  factor  to  convert  each  altitude  to  metres  (F10.3) 

H  Number  of  columns  in  the  altitude  matrix  (14) 

I  Format  of  data  in  altitude  matrix 

For  items  B,C,D,E  the  name  of  the  file  must  start  in  column  1  and  be 
followed  by  at  least  one  space. 

The  altitude  data  is  assumed  to  be  in  the  form  of  a  rectangular  matrix 
of  altitudes  derived  from  a  grid  of  points  equally  spaced  in  the  N-S  and  E-W 
directions.  The  first  row  is  taken  as  the  northernmost  and  is  read  from  west 
to  east.  If  any  zero  or  negative  value  is  encoun.ered  in  the  altitude  data 
it  is  assumed  that  the  value  is  not  valid  and  it  is  rejected  from  the  calculations. 

From  each  3x3  set  of  adjacent  altitude  points  the  following  are  calculated 
for  the  central  ooint: 

Calculated  height  of  the  quadratic  in  metres 

Aspect  angle  of  the  slope  measured  in  degrees  clockwise  from  north 
Gradient  of  maximum  slope  measured  in  degrees  below  horizontal 
Profile  convexity  measured  in  degrees/ 100m 
Plan  convexity  measured  in  degrees/lOOm 
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These  values  are  stored  for  use  by  other  programs  and  sums  of  powers  1  to 
4  ar.d  of  cross  products  are  formed.  In  those  cases  for  which  the  gradient  is 
horizontal  and  the  aspect  vertical  the  plan  convexity  is  indeterminate.  These 
cases  cannot  be  included  in  the  correlations.  Their  maximum  and  minimum 
profile  convexities  are  calculated  and  are  printed  with  a  description  of  the 
terrain.  Once  input  is  complete  correlation  coefficients  and  statistics  are 
calculated. 

Output 

(a)  Detailed  results,  output  to  'detailed  results  file'  usually  dummy  - 
output  stream  3 

Each  line  of  the  detailed  results  is  as  follows: 

Row  number,  Column  number,  Nine  altitude  values,  Calculated  altitude, 
Aspect,  Gradient,  Profile  convexity,  Plan  convexity.  For  those  points 
for  which  the  aspect  is  vertical  the  aspect  and  convexity  values  are 
assigned  impossibly  large  arbitrary  values  to  draw  attention  to  the 
points. 

(b)  Calculated  values  for  use  by  :her  programs;  these  are  output  to  the 
calculated  values  file  -  output  stream  7 

Each  line  of  calculated  values  is  as  follows: 

Row  number,  Column  number,  Actual  altitude,  Calculated  altitude, 

Aspect,  Gradient,  Profile  convexity,  Plan  convexity,  Cos  of  aspect. 

Sin  of  aspect.  These  results  are  in  a  form  suitable  for  use  by 
SPSS  and  MIDAS  and  are  used  by  the  other  display  programs. 

(c)  Class  limits  -  for  use  in  organising  display  programs  -  this  is  output 
to  the  'limits  file'  -  output  stream  9 

A  4x4  array  containing  the  maximum,  minimum,  mean  and  standard 
deviation  for  each  of  calculated  altitude,  gradient,  profile 
convexity  and  plan  convexity. 

(d)  Main  output  -  output  stream  6 

The  first  part  of  the  output  consists  of  a  list  of  the  calculated 
altitude,  minimum  and  maximum  profile  curvature  and  a  description  of 
the  type  of  point  for  all  those  points  for  which  the  aspect  was  found 
to  be  vertical.  A  summary  of  the  types  of  points  and  the  maximum 
and  minimum  altitude  for  those  points  is  printed. 
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The  second  part  of  the  output  consists  of  statistics  and  correlation 

coefficients. 

For  all  those  points  at  which  the  gradient  is  non-zero  the  mean, 
standard  deviation,  skewness,  kurtosis,  maximum  and  minimum  are  printed  for 
each  of  calculated  altitude,  gradient,  profile  convexity  and  plan  convexity. 

For  the  aspect  angle  a  vector  mean  is  printed,  and  also  a  vector  mean 
of  each  aspect  weighted  by  the  sin  of  the  gradient.  Correlation  coefficients 
for  each  pair  of  calculated  altitude,  gradient,  profile  convexity  and  plan 
convexity  are  printed. 

A  second  set  of  statistics  including  as  many  as  possible  of  the  zero 
gradient  points  is  calculated.  For  calculated  altitude  and  gradient  all 
points  are  included;  for  profile  convexity  and  plan  convexity  non-zero 
gradient  points  and  those  for  which  the  gradient  is  zero  and  the  terrain  is 
a  horizontal  plane  are  included  (for  the  'plain*  points  the  plan  convexity 
is  taken  to  be  zero) . 

Line  printer  density  map  programs  OYSM  and  OYSMA  OYSM  and  OYSMA  each 
produce  a  variable  density  line  printer  map  from  a  matrix  of  data.  Each 
data  point  is  represented  by  one  printer  position.  The  data  should  be  in 
the  form:  row  number,  column  number,  value  and  the  format  to  achieve  this 
is  part  of  the  input  data.  If  the  number  of  points  in  a  rowT  is  greater 
than  125  the  output  is  divided  into  blocks,  each  of  which  has  a  maximum  of 
125  points;  the  programs  will  accept  a  row  of  up  to  1,000  points. 

OSYM  produces  a  6-density  map. 

OYSMA  produces  an  8-density  map  from  data  representing  angles  between 
0°  and  360°. 

OYSM 

Input  Data  header  (20A4) 

Subtitle  (20A4) 

Name  of  file  containing  data 

Number  of  columns  in  matrix  of  data  (14) 

Format  of  data 

Class  boundaries  for  different  densities  (6F8.3) 


14  - 


The  name  of  the  file  containing  the  data  must  start  ir.  column  1  and  be 
followed  by  at  least  one  space.  The  data  should  be  in  form:  row  number, 
column  number,  value. 

The  class  boundaries  are  read  in  starting  with  the  lowest.  Values  lower 
than  the  first  class  boundary  and  missing  values  are  represented  as  blanks. 

Values  above  the  6tn  class  boundary  are  printed  most  densely.  The  different 
densities  are  produced  by  printing  (and  over  printing)  the  following  symbols: 

.  ,  +  ,  0  ,  0-  ,  OX  ,  OXAV 

OSYMA 

Input  Data  header  (20A4) 

Subtitle  (20A4) 

Name  of  file  containing  data 

Number  of  columns  in  matrix  of  data  (14) 

Format  of  data 

The  name  of  the  file  containing  the  data  must  start  in  column  1  and  be 
followed  by  at  least  one  space.  The  data  should  be  in  the  form:  row  number, 
column  number,  value.  The  values  are  assumed  to  lie  between  0  and  360.  Any 

values  below  0  or  greater  than  360  are  treated  as  if  they  were  0  or  360. 

Missing  values  are  represented  by  blanks.  The  values  V  are  assumed  to  be  measured 
clockwise  from  N.  The  most  northerly  are  represented  by  the  most  dense  printing 
and  the  southerly  values  by  the  least  dense.  East  and  West  are  of  similar 
density  but  use  different  characters  as  shown  below  , 


North  V  <22.5  and  >337.5  OXAV 

22.5  <  V  <  67.5  X+ 

East  67.5  <  V  <112.5  X 

112.5  <  V  <157.5  + 

South  157.5  <  V  <202.5 

2C2.5  <  V  <247.5  0 

West  247.5  <  V  <292.5  0- 

292.5  <  V  <337.5  OX 


Slope  plotting  program  OYPL  The  slope  plotting  program  OYPL  is  used  with 
the  MTS  plotter  routines  in  *PL0YSYS  to  produce  a  map  in  which  the  magnitude 
and  direction  of  the  gradient  at  each  grid  point  are  represented  by  an  arrow. 
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The  map  is  29"  wide  and  the  distance  between  pic  teed  grid  points  depends  upon 
the  number  of  columns  in  the  matrix  of  data.  The  shaft  of  the  arrow  is 
centred  on  the  grid  point,  and  is  in  the  direction  of  maximum  gradient.  The 
length  of  the  arrow  depends  upon  the  magnitude  of  the  gradient  and  can  be 
zero  or  from  1/5  to  5/5  of  the  distance  between  grid  points.  Five  class 
boundaries  are  read  in  and  if  the  magnitude  of  the  gradient  is  less  than  the 
first  class  boundary  only  an  arrow  head  is  produced;  if  the  gradient  lies 
between  the  first  and  second  intervals  the  arrow  is  1/5  of  the  grid  length, 
etc. 

Input  Data  header  (20A4) 

Name  of  file  containing  calculated  values 
Number  of  columns  in  matrix  (14) 

Class  intervals  (5F6.3) 

The  name  of  the  file  containing  the  calculated  values  must  start  in 
column  1  and  be  followed  by  at  least  one  space.  The  program  reads  row 
number,  column  number,  aspect  and  gradient  from  the  calculated  values  file 
produced  by  the  main  program. 

Class  limits  program  OMYL  The  class  limits  program  is  an  organisational 
device  to  ensure  that,  for  a  wide  variety  of  altitude  matrices,  acceptable 
graphs  and  maps  are  produced  automatically.  The  program  produces  scales 
and  limits  for  use  by  other  programs  and  writes  these  into  appropriate 
positions  in  a  master  file  of  run  commands  and  data.  Output  stream  7  is 
initially  set  to  the  master  file  (GEOD) . 

Input  Name  of  file  containing  limits 

The  limits  file  produced  by  the  main  program  holds  a  4x4  array 
containing  the  maximum,  minimum,  mean  and  standard  deviation  of  each  of  the 
following:  calculated  altitude,  gradient,  profile  convexity  and  plan 
convexity. 

The  MIDAS  histogram  program  requires  a  minimum,  maximum  and  class 
interval  for  each  histogram.  For  calculated  altitude  which  is  the  most 
likely  to  vary  for  different  matrices,  the  class  interval  is  chosen  to  give 
between  60  and  120  classes  and  the  maximum  and  minimum  are  rounded  up  and 
do™  respectively  to  suitable  values.  The  minimum,  maximum  and  class 
interval  are  then  written  to  a  specific  line  in  the  master  file  in  a  .format 
suitable  for  the  histogram  program. 
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For  the  line  printer  density  maps  variable  limits  are  used  for  altitude 
and  gradient.  The  class  boundaries  used  are  (where  m  is  mean  and  sd  standard 
deviation):  m  -  1.2sd,  m  -  0.6sd,  m,  m.  v  0.6sd,  m  +  1.2sdjand  m  and  sd  are 
rounded  to  give  suitable  class  boundaries. 

Output  The  output  consists  of  lines  of  data  written  to  appropriate 
lines  of  the  master  file  for  use  by  other  programs. 

Master  file  GEOD  GEOD  is  the  master  file  containing  all  the  run  commands 
and  data  necessary  to  run  the  terrain  analysis  programs.  The  run  commands 
start  at  line  -1000. 

Lines  1  to  9  of  the  file  contain  the  file  names  and  data  pertaining  to 
a  specific  run  and  are  used  from  there  by  several  programs.  Lines  1  to  9 
contain  the  following: 

1  Data  header  (20A4) 

2  Name  of  file  containing  altitude  data 

3  Name  of  file  for  detailed  results  -  usually  *DUMMY* 

4  Name  of  file  for  calculated  values 

5  Name  of  file  for  limits 

6  Grid  size  in  metres  (F6.3) 

7  Scaling  factor  to  convert  altitude  data  to  metres  (FI 0.3) 

8  Number  of  columns  in  altitude  matrix  (l4) 

9  Format  for  altitude  data 
END  OF  FILE 

REFERENCES 

(1)  Peucker  T.K.  and  Douglas  D,H0  (1975),  ’Detection  of  surface- specific 
points  by  local  parallel  processing  of  discrete  terrain  elevation  data’, 
Computer  Graphics  and  Image  Processing  (1975),  4,  (375-387). 

(2)  Grender  G.Co  (1976),  ’T0P0  III  :  A  FORTRAN  program  for  terrain  analysis’, 
Computers  and  Geosciences  (1976),  2,  (195-209), 

(3)  Evans  I.S.  (1977),  'The  selection  of  class  intervals’,  Trans.  Inst,  of 
British  Geographers  (1977)  N.S.2,  1,(98-124). 


Figure  1  Specimen  of  input  to  GEOD. 


Figure  2  Statistics  produced  by  main  program  0MY8. 
Data  for  Torridon,  N.W.  Scotland. 


Figure  3  Part  of  output  of  main  program  0MY8, 
listing  points  with  zero  gradient. 
Data  for  Torridon,  N.W.  Scotland. 


Figure  4  Output  of  line  printer  density  map 
program  OYSM.  Calculated  height 
for  Torridon,  N.W.  Scotland. 
Hissing  values  have  zero  gradient. 


Figure  5  Output  of  line  printer  density  map 
program  OYSM.  Calculated  height 
for  Cache,  Oklahoma.  Missing 
values  have  zero  gradient. 


Figure  6  Output  of  line  printer  density  map 
program  OYSM.  Maximum  slope  for 
Torridon,  N.W.  Scotland. 


Igure  7  Output  of  line  printer  density  map 
program  OYSMA,  giving  aspect. 

Data  for  Ferro,  Calabria,  S.  Italy. 


Figure  8  Output  of  plotter  program  OYPL. 

Data  for  Nuppr,  N.W.  Iceland. 


Figure  9  Histogram  produced  by  MIDAS,  showing 
profile  convexity  for  Torridon,  N.W. 
Scotland. 


igure  10  Scattergram  produced  by  SPSS,  showing 
calculated  height  against  gradient  of 
max.  slope  for  Torridon,  N.W.  Scotland. 


Figure  i  Specimen  of  input  to  GEOD 


■#END 

# 


.1  TOPS  I  WON  SCOTLAND 

2  GUN  :l.  t  T ORR  ( 4  ?  303 ) 

3  HtDUMMY# 

A  TORSE V 

5  INTT 

6  100.00 

7  1,0000 

8  100 

9  i 20 FA . 0 ) 

OF  FILE 


IHIS  PAGE  IS  BEST  QUALITY  PEAQTIAABLE 
•I5R0M  COPY  FISSIISHED  TO  DD.C  _ _ 


£1^ure-l  Statistics  produced  by  main  program  01IY8 
.Data  for  Torridon,  N.W.  Scotland. 
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Figure  3  Part  of  output  of  main  program  0Mf8, 
listing  points  with  zero  gradient, 
Pata  for  Torridon,  N.W.  Scotland. 
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Figur.  4  Output  of  line  printer  density  tup 
progran  Oi'SH.  Calculated  height 
for  torridon,  H,U,  Scotland. 
Missing  values  have  aero  gradient. 
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Figure  S  Output  of  line  printer  density  mop 
program  OYSM.  Calculated  height 
for  Cache,  Oklahoma.  Kissing 
values  have  aero  gradient. 
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riture  6  Output  of  line  printer  density  oap 
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M.U.  Sratland. 


$0 


rivur*  »  Output  mt  »Utt«r  fio|M*  CWt. 

feu  f#r  K.V.  Itplmi. 


NUPUR  NW  ICELAND  100H  16  AUG  1977 
SLOPE  VALUES  1.000  2.000 


10.000 


15.000 


xV»\NN*v» 

NSSSV.- 

SSNNN-*- 


a  pi 

I  7  ’  7  f  VSiSS*.***  <  W  /  <  “  V\N 
*  ♦  4  4  4  4  SVVS 

T..  14!  I1.  ..iiujj 


. r  V  \  \  \  \  J  1  1 

■!  t  i>!<: : !  iaSssssSJSSS^'l^^^v?^ 


I  ?/'*'*■*  *  * 
7/S"*** 
/"*•"*  * 


— • ' '  m  l 

i  ! !  I^wHJJvSvSS; ; 

mk^sdMm  |  ji^awp 

♦ <U \) )K. : : 1 i yS&vg?  f f  Ucs™ 


- - - - - -  T  ■\"S'WV>  \  N  N  \  \  ' 

'  NSvIv^>-  -  *  '  VKVN\WV*sV*s  \  //// 

■  •  i  .  ^  WSSSSVV*^1^  7  i  WWWSSV^  7 

^VSSSNNNNV^^J  |  V\\\I$is^U\SY^ 


>// /« » •<#* #  #  m  i  v i 

✓,*/44  4  \SV\l  j>V»4/ 

a3f«i«i  i 


4 \SSSSNS7{ 7 

/SSSSNVf  I 


VsY/l \\\\\*<**'*  •'*-*-*  •  t 

//^//{sSU/7.  «  •  *7W.  •  •  • 

'ssxc\ sssi^-4 :  sn!/££»  i  *  - 

Insww'  •*** 

CCwt\N\/Iv/  IssiWi  **** 

L'.XCZZ&i  rauysAi  <«<  \  j 


mm  m 


(fcsssi WM  :  U 

*  >  /  /  ^*TT  In  I  7—  *  vrr»  t  »WV<-rrV*^  1  1  J  X  J 

J  -w.  •  ^///-^//  J  /✓W'V'."  .  <»<*(•  !'*^jif  i'Vj  ,  .  .- 

V  -V>>/  ///*A*->  •  *  /VVV'V'VV” *  ^  i  i  ?  f  i 


1SSSSN- 

SN»»: 


::S"i£K5j|  If® 
!{!•#{£: 


.  /  /-•*-—. 

►//«**' 

rvv*-!/ 


WZ£Z%%%  te ;•:::! SiS n ! III! 

fell  H!1H1«^MI 

! II! . : 1 u  4 u 4 n  ?7  7#/vv<» 

!?;?  :!  .. •  -'•vy'S/i  1 1 l J 4 1  'tWSsZS/ji  *  *  * 


ff>;^ 


<#>77. 

M  M^i 


:  i  i  H  i  J ! } !  U  j :  i  i :  ss 

: :  5a« : ; ; ;  5  5  5  1  in  ;  t : 


I  fc«!  5  ^  !  7 


i  / / 4  <** *  *  *  * 


R2»^U 

?sv.NS\y 


I  ?A.V/<  'NW 
J  ■»“*-*- J~r- 

1  /rw^— ..  •/A*/*-' 
4  .-^.-4  .  /A*^« 


-*.\J  j 

— NJJ 


4  t  ‘ 

r  l  »-«^v'<  4 1  v 

^llltUllM  Sljjjiiim^ 

:  .  :::HSB332sii5 iUU ^ 
Scg22:i : ; . ; : 

»<>♦♦*»♦•«*»  M>r '  *  >  i 


*  *  '■ /  < '. '. '.s..  ...j 
HI  1  i  \  -iSSftf/S-* 


7i<<Wi»555.r  U' 


ill}! 


\\  1  /// \SNN»S<V»'»»  »  ••*• 

U  /  // i  {nnnnv*  »•*“ 

U///1  «— - 


*«  ^iu «  i 


♦  ***•*»  »  <>V7. 4  s\\  J ~!*vCC^ 

tin....  *-*7\  »  »  •••'.'•'•Si  1  V  W'C' 


f; 

a 

o 

o 

o 

o 

u 

o 

o 

o 

O 

O 

■f 

o 

u 

o 

o 

c 

o 

o 

o 

o 

o 

Q 

‘•i 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

m 

* 

m 

• 

«<» 

m 

« 

r» 

«» 

ru 

m 

•» 

cr 

c 

o 

H 

«4 

IM 

IX 

m 

•> 

c 

>0 

»r 

in 

fst 

m* 

•  ♦ 


a  » 

rv 

» 

IM 

l  a 

c 

a-  1 

«l 

Ml  | 

♦ 

«v* 

c 

•> 

«OI 

« 

♦ 

IV 

« 

«  aa 

0 

*  «  ft 

ft 

« 

ft 

m  1 

ft  « 

a  1 

• 

W  | 

V# 

ft  ♦  ♦ 

ft  « 

a 

• 

* 

<M  4 

*  ft  ft  ft  ft 

i-  1 

■4  ft  ft 

r» 

4  i 

#  % 

ft  »  ft 

ft  ft 

V 

«*i 

SI 

♦ 

♦  *  * 

4*  * 

ft  *  ft 

ft 

Mft  ft 

c  run  *  * 

rv* 

an 

» 

*  Mb  * 

M  ft  ft  ft 

rv* 

V 

a  4  i 

«  ft  * 

II  #  CMfvJ 

ft  rvjrvjrvi 

♦ 

IU 

C  } 

mm  r  r 

•  ft  ♦  ft  ft 

ft  ft 

u 

i 

» 

ft  fvift  ft  rsm*  «  N 

IN* 

23  S 

u  u 
o  C  O 
XI  •  u 
•  *H  (A 
■w 

•  «  • 
tfi  u  tt 
ft>  (A  • 
tu  SB 
</)  u 

l»  C  B 

B'j  S 
•0&5 
S  *  fc 
£  a  £ 
IV 

B  >H  h 
A  II  O 

g** 

«•  T*  •» 

B  II  A 

lf«5 

«l  w  « 

as  . 


i»  #ft 

ft  Aft 

IV 

*ft 

*M<V«  | 

ft 

M 

U  | 

» 

*1  M 

▲  1 

N 

♦ 

M  ft  - 

c  1 

ro  ■ 

1 

*  ft 

• 

* 

ft 

N«i  |4 

11 

I'QUI  ft 

ft 

# 

4 

IW.MJ 

a  »  ft 

ft* 

N-  ft 

ft 

Vft*4  ft 

ft 

•N 

«4*W  ft 

ft 

ft 

V 

II  1 

ft 

« 

2» 

•  I 

• 

•» 

v 

Wa 

**■  •• 

ft 

•ft  V  ft 

m 

W  1 

IM  ft 

1 

ft 

ft 

afv.  ft 

*  ♦ 

C 

4*  ♦ 

V* 

Q 

ft 

M 

r*ft  « 

vildft 

c. 

U*W4  1 

a 

•i 

uw  aft 

W*M»t  | 

w 

*  ft 

•  u\<ri>ia 4tna « »  w  «  o 

•  a « ivrviv<*iin*vra<v  pwaa  impi*  p'  M" 

•  a  a •  anipimivripim# •  iva  i*m  X  a  • 

fna  #«vo  nirvuipiin#  tnPa4PaivtV44  S  g  . 

«  •  «»|P«m«  P.*#m4 P>a  4  Pa  a  ana 

•  H  IPIVIVNIVIV*  vCM4»  «M«  «<  «  4  4 

«  *  «  <  Pi  » tV4PlM&4IMn4  V  4N4'  a*  mob 

a  4  in  *4414  « rvar>iu\4iMiM*M4aiir  ««« 

r  Mr  •nyaiiti  iv®4iMipi"tni*i4r-i*i*  # 

a  n.ivrva  *<*‘1!  p.iviMvin4PJ4P>PiP4PHi  in  p*  o 

IV  4.  pa«IVPip.|(VtCPiP>4Pipa(V.NinlV4<iJ«  4*  ¥  — 

iv  an.  nia  b  a  *  nvriaimojv  na  p.  v#  a  «  a  a 

«V4ls*4#4P*Pa#  VlVIV4*IV4»ipiPiPi#  «ui  JJ 

t  •  a*ll  #1*1441114  44P4PI|ftPI4NPaPllVa  4  3 

a««  IP  **i4Pia  |I|IVMI>I«  4<apwMffl4i  iv# a  M 

a  Pa  pnv4iv<vM»0  B>P'fviv»iv4y  i>iivra4«  *  •>* 

«w»«ivi>,tvt  «aa  w  rvi«ia«*»v4iv»4#pa  iviva  to 

*  iv#t>i«t'a«iv4i*iP'faivpt«.«ift4iit<«<Papi  ♦P'l" 

W  •«  IV#  4  IVHI  IV«  lVfVP,IM44«V*P,PI*  a  p* 

*  a  »  p.#a«4  4iv«vpnvpiiv*4a'«VP»*P'lva*  a 
a  aiviv  pi«  a«viviv«  Pia<viv4m«r>«  aiv«iv«« 

«v  at  piauv#  4*v# a  ivpipimv'*  oiaapuvraiva  r> 

a  iv  a  *  ivp>  a  *  m  a  vmvii  iM«  iurva  iv#  ogtvtnr  a  piiv«  p 
a  p  iv  iv«*  iv#pu«i»  <v<vra«  pi#p>4#Pi4a>a  miv  • 

*.  a  a  a  *  k4ttPta4P'iM«tvoi(V4<V4  4**ra<*w\  iv  w 

A  a  P*Pa«  «  rvy  w  Ivray  pimpi  y  MPif<pinP.#4  X 

x  *  a  a«ra  -»  im  a  v  v  mivojin  a  x  4fvra4  jnani* 

iv  iv  irnvx  *a>i*iP'4ivMPipiw<4iio4P*»v*aiv  a  * 

#  •  »» a  1*101  1*1*  «  4P>tt*vpi*  VM4UIVPI44PHV*  P*»  X 

,  «  ivp*  *  pa  pi#  «  4»*<r*'4rifl*>«  inn  a  a  a  * 

.*  a  «  .£tv  ivraNvir>i*i«PaivN*risu>pa4<rinii\PiPii*)«  a  «v*. 

pa*  i-i*a*  4# in*  «  »Mi  «*o>*na,'.*44piivx  *  a  « 

♦  a  iva  « IV4PIPIN# «pi  ror-4)4*uin4Pi«*i<ip>  m  w  aiv« 

*  a  p*  4  a  a  w«  *  4v*v#a  iv#  #  a>*gP>«i«*Piui*ii<v(v. 

a  w  pa  m4iv  *  p>mmiPrPivin<gift<gtm44V  44  *v  p>  «  y  n 

a  a*  •  a  iv  a*  iv  a  iV4**uv  a  »**N444Wm4  a  oiptaiva  4#  »v  v 

a«  iv ,  a  iv  i*ia  ivp-piapnvivivp*iwrvaipiMU'MBiviv#«pia  a  *v 

«  pi  »  « pia  v  ivp'pi#  a  ma  iv«ivirPHr<fM*'t*4rvir4p>4  miviv 
a  aa  **  n«40>  p>ain<gpicv*v4iv,*4*oino>  o»«y.<gn  *a  <v 

ro  aiv  ni  # « n>inivpiv  uik  p>4i*>ini~u|iri«(r4iviiMvpipi4  a 
pi  a*  rviv  mi'nvivuia « u>u\4  voivir4ioiKP*oaivv<Pia t*i*v*» a  iv 
«*#  a  rviv*  a  ivnunw  atptv4*pp<(v*>4r»g'»o>na4aipi«vP'4  #iv  ca  rv 

a  a  an  a  <v  iv  p-  a  rvni4  a  iniMpj4aua'ip«V'if  va*o*P-p  rv  ir«v  ,  > 

•  *oiiv  4<va  pipipi  ivpiPir*>iriM*i9><*i«Prvv*4V>Pi  iv  eg. 

p  I  «  a  #  a  a  iva  4*  4#  iviva*ai44HP»*440»«ivaa  o»#p>i*i  a  a» 

P'  rvrva  av»  a  <t  4‘V(viP4<gg,trv*wipiPii,«>4«VMi*  a  ivp*  a  t 

a  #*  P.IV#  4MIV0'4  4  444PIP>I»V‘V'V  <44V'l/'V*0'P»VHVPIf  *al*.|\ 

a*a  rvaaiva  n«voiU'Pi40it.  M»4*vwu,tfW,<*»iPiPV,4>P>P*Pi4Pi#  iv«  tv 

#a  iva  •••iii4iwii*ipii**ojvi«>aiV'V'«rv‘ipw'V*oiV  i**i*iviivrva  a  iv 
i|  a  a  iv  a  ««4*  ap»4<*iui»lPr»i«ipu'li*lM*’V,W'W,V,**l“P*i”apia  a  u*  a 

pi  iva  v  tv'*l4*>wr»4a>avg*f*r»w<v  v'i»’*v*»*r,4«*,'i»ia  iv 

a  a  4ini*i4P,,vonuv'v'4>Cv'W'**iiMi  v«'vrv  a'aip'4*”*i"4  » *  * 

*  * » a  iviv  #  i**iv  <**lvv>v'~'*V|vv'|pv*U'UiU'V'V'V  *#V'l/*4i**pi*  a 

a  'i  a  ^ivir-o't-4in4h-4o>oiirw  av'CMi'u-Q'iniMMAmoi  piipiv  n 

a  a  y  «v*  v\>caMV4ing'g  <av V'tf,m<B«'4frV'V  V'v,oiv  *v»'iM  a  a  p 
a  «  *  oi4inB''*voiv  g  u  wv*g  vpg  v  g'V'g'g-g  g',v*».*vi*itv  a  a 

a*  a  a  ni4ojrvaptojp*4o}'00'g'g  u>g'v>vg'i/>g’VV,g  #4*v# 

im  «  r  «  •  npip*4V  <og*g-«i«'w,g*g,g  w  g'g*w*g*g,g'4v*ivv 

«  aa*  a  a  iv  a  o*o;<o«>*(ig'g*<ti<*>w>g'g'fpg*v*w'g  g'g*w  4  iviv  iv 

*  •  i  i  PiPUvMlt  *  iPinw«*«rg'i<*P*4g*iripg  gg'g  g*vP-Pia  4  v 

a  •  »  a  ♦  m»  *Nr'aiM*iirtf  oiiBV’ga'ii'VvvB  vg'«  v<t4Mi 

>*ia  a*  iviv*  i*»*ia  vuiip-l»i*i'*1,  a  *gg'g,g‘g'g,v>g’g*g*g,w'iv4  aiva 
a  mi  pam4  a  a  y  g  ipg'g'iPM'gg’g'V'g  g*g'«.  ivin  * 

a  a-  ivivipivpipi  « ivu>w*irg'g'v>g*ipg*g  ipi»iv  » . 

a  »  <*  iv  h  ivo4iw*i‘iP*4i*ig  V'ipg‘g,g,g’«,g,g  g*w  a  •* 

a  ivcva  a<vivip*g«w  #  4U*o*g*g>tpg  g*g*g  g  aia  aiva  4 
*pp  «  «*ia*g'vgw*gi4wgg'v  in'*  a 

a  a  *  a  #4i*i*  a  a  m®444Pmg,iMpg  44*  «v 

a  »r>piiv  *•  4Hg»4V'««Pipw*u<ip 

a  a  aivpnviva  4  ivivpi  a 


V* 

S 

U*u 

ft. 

ftift 

MW#M»«MN4NNMWrwNNMr»<IMMHrM(i|NMtMM 

-P+++  ♦  •’*►4*4 

^♦WNN 

»-♦  ftO 

« 

u 

W 

ft 

o 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

W4ft 

u 

W 

W* 

ft 

ft 

vft 

ft 

ft 

ft 

ft 

ft. 

• 

• 

ft 

ft 

ft 

• 

ft 

ft 

ft 

ft 

ft 

w* 

4* 

ft* 

MU 

IM 

u> 

Ol 

a* 

1*1 

a> 

•*l 

MW 

0* 

w 

ft* 

«* 

ftp 

ft* 

ft* 

w 

O 

*4 

IV 

r* 

04 

« 

u 

ft* 

Uft* 

ft* 

ft* 

ft* 

>o 

m 

1*1 

IV 

•4 

1'HIS  PAG£  Xi  QUALITY  PRAC3 


